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Single- and two-particle excitation spectra of the one-dimensional, half-filled Holstein-Hubbard 
model are calculated using the continuous-time quantum Monte Carlo method. In the metallic 
phase, the results are consistent with a Luther-Emery liquid that has gapped spin and single- 
particle excitations but a gap less charge mode. However, given the initially exponential dependence 
of the spin gap on the backscattering matrix element, the numerical excitation spectra appear 
gap less in the weak-coupling regime, and therefore resemble those of a Luttinger liquid. The Mott 
phase has the expected charge gap and gapless spin excitations. The Peierls state shows a charge, 
spin and single-particle gap, a soft phonon mode, backfolded shadow bands and soliton excitations. 
Arguments and numerical evidence for the existence of a nonzero spin gap throughout the metallic 
phase are provided in terms of equal-time spin and charge correlation functions. 

PACS numbers: 71.10.Hf, 71.10.Pm, 71.45. Lr, 71.30.+h, 71.38.-k 



I. INTRODUCTION 

The interplay of electron-electron interaction and the 
coupling of electrons to the lattice is a fascinating topic 
in condensed matter physics, with direct experimen- 
tal relevance in compounds such as transition metal 
dichalcogenides.^ Theoretically, this problem is partic- 
ularly difficult in the case of non- Fermi liquids, e.g., 
Luttinger liquids realized in quasi-one-dimensional com- 
pounds. The difficulty is associated with the competi- 
tion between the instantaneous Coulomb repulsion and 
the retarded, phonon-mediated attractive interaction. 

The Holstein-Hubbard model is one of the most 
frequently studied models in this context. It com- 
bines a local electron-electron repulsion (with strength 
[/), a nearest-neighbor hopping and a local density- 
displacement coupling between electrons and phonons. 
Even though it neglects anharmonicity effects and cor- 
rections due to incomplete screening,^ it is believed to 
capture the key physics. In one dimension, where power- 
ful numerical and analytical methods are available, exist- 
ing work has mostly focused on the bipolaron problem of 
two electrons, see Refs. 3 and 4 for an overview, and the 
half- filled band.l^HI^ Some results fo r intermediate band 
fillings are also available .-^'^^^^^^'^ For vanishing Hub- 
bard repulsion, the Holstein-Hubbard model reduces to 
Holstein's molecular crystal model^^ which has been ex- 
tensively studied at half filling, see Refs. 2 22 - 25 and ref- 
erences therein. The Holstein-Hubbard model may also 
be investigated in the framework of dynamical mean-field 
theory,^^ although non-Fermi liquid physics and spa- 
tial correlations are not captured. 

Interaction effects are dominant for a half-filled one- 
dimensional system. In the Hubbard model, a repul- 
sive electron-electron interaction leads to a Mott insu- 
lator with dominant spin-density-wave correlations but 
no long-range order (a detailed definition of the various 
phases considered in the following is given in Sec.[lV| see 
also Table |l|. If the interaction is attractive, arising from 
electron-phonon coupling in the limit of high phonon fre- 



quencies, backscattering opens a spin gap , b ut the sys- 
tem remains metallic. Hirsch and FradkiiP^ argue that 
the U{1) charge symmetry is broken at finite phonon fre- 
quencies, so that long-range order can exist at T = as 
a result of the Peierls instability.^^ Alternatively, such a 
transition to a state with two electrons at every other 
site can result from a nearest-neighbor electron-electron 
repulsion that emerges from electron-phonon interaction. 
In the static limit of classical phonons, exactly solvable 
in mean-field theory, an insulating Peierls state exists 
for any nonzero electron-phonon coupling. To fully un- 
derstand the role of quantum lattice fluctuations (a fi- 
nite phonon frequency, as relevant for experiments) has 
turned out to be a complex problem that still poses a 
number of open questions. 

Whereas early numerical results for the half-filled 
Holstein model suggested the existence of an insu- 
lating Peierls state for any nonze ro coupling, in 
agreement with stro ng-co uplin^p^^ and renormaliza- 
tion group arguments p * ^^ * subsequent and more accurate 
treatments strongly suggest that lattice fluctuations de- 
stroy the charge-ordered state below a critical coupling 
strength. Turning to the Holstein-Hubbard model, a 
metallic phase of the Holstein model should survive for 
small enough val ues of the electron repulsion, as pre- 
dicted previously.'^^ Early numerical studies focused 
on large values of U where the system is only metal- 
lic at the single point of the Mott insulator to Peierls 
insulator transition. Evidence for a metallic phase, 
whose extent stro ngly d epends on the phonon frequency, 
was reported later .'^"^^^ Recently, work aiming at char- 
acterizing this intermediate phase has initiated a de- 
bate about the p ossibility of dominant superconducting 
correlations the existence of a spin gap,^the valid- 
ity of Luttinger liquid theory for retarded interactions,'^ 
and the existence of a metallic phase."^^^^^^ 

The metallic phase is of particular interest in relation 
to superconductivity in quasi-one-dimensional or higher 
dimensional systems. Even though the continuous U{1) 
gauge symmetry cannot be broken in one dimension, 
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dominant superconducting correlations imply a tendency 
toward superconductivity in higher-dimensional settings. 
Initial numerical results for the metallic phase in the 
Holstein-Hubbar d model indicated a Luttinger liquid pa- 
rameter Kp > i| 8 | 9 | ii | i2 ] suggests dominant pair- 
ing correlations. However, direct calculations of correla- 
tors instead reveal that charge-density-wave correlations 
dominate. I^^inil This contradiction has been attributed to 
finite-size effect^ and logarithmic corrections^ due to 
the retarded nature of the phonon-mediated interaction. 
Dominant pair correlations together with Kp > 1 have 
been reported at quarter filling, and degenerate pair- 
ing and charge correlations exist in the half-filled lattice 
Frohlich model. ^ 

Another controversial issue is the existence of a spin 
gap in the metallic phase, as a result of binding of elec- 
tron pairs into singlet bipolarons. This would imply that 
the low-energy theory is that of a Luther-Emery liquid (a 
Luttinger liquid with gapped s pin ex citations but gapless 
charge excitations, see Table Ib,'^^ which has also been 
suggested as the appropriate description of the normal 
state of Peierls compounds (i.e., for T > Tc)^ In the 
nonadiabatic limit, the mapping to the attractive Hub- 
bard model implies a nonzero spin gap. For small in- 
teractions, the gap scales as^^ A^r ^ e~^^/^ where 
the effective backscattering matrix element, can either 
be related to the attractive Hubbard interaction medi- 
ated by the phonons in the Holstein model, or to the 
net interaction resulting from the interplay of electron- 
phonon and electron-electron interaction in the Holstein- 
Hubbard model. In the static mean-field limit, a spin gap 
exists for any nonzero coupling. For finite and especially 
for low phonon frequencies, there are significant retarda- 
tion effects. Similar to the destruction of the Peierls state 
by lattice fiuctua tions due to a reduction of umklapp 
scatteringj ^^ l ^-^ l ^^ the spin gap may be destroyed for small 
enough interactions by a renormalization of backscat- 
tering under the renormalization group fiow. Another 
possible mechanism is the dissociation of bipolarons (a 
bound state always exists for exactly two electrons, but 
this bipolaron is spatially extended for weak coupling) 
due to mutual overlap, similar to polaron dissociation at 
finite band filling. On the other hand, in the low- 
energy limit, any finite phonon frequency may naively be 
argued to be irrelevant, leading back to the spin-gapped 
attractive Hubbard model. There is also evidence from 
renormalization group calculations that retardation can 
enhance backscattering.^^ 

A gap that opens exponentially is of course very dif- 
ficult to detect by numerical methods. Several previous 
works explicitly state th at the metallic phase presum- 
ably has a finite spin gap.^'^^^ A crossover from gapless 
to gapped spin excitations inside the metallic region has 
also been proposed. Recently, a spin gap has been ob- 
served for all phonon frequencies in the quarter-filled Hol- 
stein model at intermediate electron-phonon coupling.'^ 

Here, the continuous-time quantum Monte Carlo 
(CTQMC) method is used to calculate the single-particle 



spectral function as well as the dynamical charge and 
spin structure factors in all three phases of the Holstein- 
Hubba rd mode l. Previous work focused on single-particle 
spectraE^^l^^'^ and the optical conductivity.^^^ Addi- 
tionally, arguments and numerical evidence for the ex- 
istence of a spin gap in the entire metallic phase are pre- 
sented. The paper is organized as follows: After briefiy 
introducing the model and the method in Sees. [ll| and |III[ 
the results are presented in Sec. \LV\ followed by the con- 
clusions in Sec. El 

II. MODEL 

In one dimension, the Holstein-Hubbard Hamiltonian 
can be written as 

^ = E + H.c.) + [/ ^ n^tni; (1) 

ia i 

+E(i+f^)-»x:*(n,-i), 

The first and second terms constitute the Hubbard 
model, describing electrons with nearest-neighbor hop- 
ping t and local repulsion U. The third and fourth terms 
correspond to the lattice degrees of freedom and the 
electron-phonon coupling. The phonons are described as 
harmonic oscillators with frequency ujq = ^/K/M^ dis- 
placement Qi and momentum P^, and the coupling is of 
the density-displacement type proposed by Holsteinf^ g 
is the coupling strength. The usual definitions fiia- = 
4^c-^, Hi = fiia, and n = (fii) are used. 

Using the path-integral representation, an effective, 
frequency-dependent electron-electron interaction 

U{u) = U^XW^^ (2) 

can be derived, where A = g'^/{KW) is a dimensionless 
electron-phonon coupling constant and W = 4t is the 
free bandwidth. In the nonadiabatic limit (high phonon 
frequency, Uq oo) or at low energies {u 0), this 
interaction reduces to an instantaneous attractive or re- 
pulsive Hubbard interaction Uoo = U — \W . In the fol- 
lowing, the hopping t is taken as the unit of energy, and 
the lattice constant and h are set equal to one. 

III. METHOD 

The model ^ can be studied with the CTQMC 
method in the interaction-expansion formulation.^ This 
method is free of Trotter errors, and has been success- 
fully applied to electron-phonon lattice models.'^'^^ESEIl 
The bosonic degrees of freedom are integrated out exactly 
(without a cutoff for the bosonic Hilbert space), and the 
resulting fermionic model with retarded electron-electron 
interactions is simulated. Methodological details can be 



found in previous publication^^I^^E^I^IIini ^nd a review!^ 
The numerical effort scales with the average expansion 
order, making large system sizes accessible at weak cou- 
pling. The method gives exact results with statistical 
errors for imaginary- time correlation functions. 

The results below have been obtained at low but finite 
temperatures 20 < pt < 162; the inverse temperature is 
specified in each figure caption. When comparing differ- 
ent system sizes, the ratio I3t/L has been kept constant. 
The system size was L = 30 or L = 50 for excitation 
spectra, and L = 30-162 for static correlation functions. 
Each QMC run was started with a zero-vertex configu- 
ration, and equilibration was carried out until the accep- 
tance rates for addition and removal of a single vertex 
have become equal to within one percent. Measurements 
were collected from independent, parallel runs (typically 
31) for each set of parameters. For the most demanding 
parameters considered (Fig. [t]), about 500 bins of 150 
sweeps each were recorded (with measurements taken at 
the end of each sweep). A sweep corresponded to 1000 
proposed Monte Carlo updates (addition or removal of 
a single vertex, or 32 attempted flips of auxiliary Ising 
spin^^). Error analysis included the usual binning and 
jackknife procedures to avoid underestimation of statis- 
tical errors due to autocorrelations. Error bars shown 
indicate the standard error. 

Spectral functions were obtained by applying a 
stochastic maximum-entropy methoci^to the imaginary- 
time Green function data and their covariance matrix. 
Convergence of statistical errors was achieved over at 
least three orders of magnitude. A flat default model 
was used, and the renormalization parameter was cho- 
sen to (approximately) achieve = where L^- is the 
number of points on the imaginary time axis. 

The static correlation functions of interest are 

Sp{r) = {{fir - n)(no - n)) , (3) 
5f(r) = (^,^^o^), 

P(r) = (At Ao) (A. = 4^4^), 

measuring charge, spin, and s-wave pairing correlations, 
as well as their Fourier transforms. In combination with 
the maximum entropy method,!^ excitation spectra can 
be calculated, including the dynamical charge structure 
factor {Aji = Ej — Ei)^ 

Sp{q,uj) = 1 ^ I (i| \j) fe-^^^5iA,i - cv) , (4) 




the dynamical spin structure factor, 

SAq,^) = I E I 1^') l'e-''^^^(A,, - u) , (6) 

ij 
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FIG. 1: (Color online) Schematic ground-state phase di- 
agram of the half- filled, one-dimensional Holstein-Hubbard 
model with a finit e phonon f requency, as suggested by nu- 
merical simulations.'^™'^^ For U/W > A, the system is a 
Mott insulator with dominant, power-law SDW correlations 
and a nonzero charge gap Ap (the notation COSl is explained 
in the text and in Table |l|. For A ^ U/W, the ground state 
is a Peierls insulator with long-range 2A;f charge correlations 
and nonzero charge and spin gaps A^, A^-. For small U/W 
and A (the scale depends on the phonon frequency), a metal- 
lic phase with gapless charge excitations exists for A > U/W 
{Uoo < 0). Here it is argued that the entire metallic phase 
has a spin gap. 



and the single-particle spectral function 

ija 

(7) 

Here \i) denotes an eigenstate with energy Ei. 



IV. RESULTS 

Figure [l] shows a schematic phase diagram of the 
half-filled H olstein-H ubbard model. Exact numerical 
investigation j^ l ^ l -^^ ^ suggest the existence of three dif- 
ferent phases at finite phonon frequencies. These phases 
can be characterized and distinguished by the presence 
or absence of an excitation gap for single-particle (Agp), 
charge (Ap) and/or spin excitations (A^r) with zero mo- 
mentum, see Table [ij Following Balents and Fisher the 
notation CxSy is used, where x (y) is the number of gap- 
less charge (spin) modes (with x,y = 0,1 for a strictly 
one-dimensional system). Since a two-particle gap for 
either spin or charge excitations also implies a nonzero 
single-particle gap Agp, knowledge of Ap and A^r is suf- 
ficient. Finally, because low-energy charge transport is 
determined by long- wavelength density fluctuations, the 
charge gap Ap further distinguishes metallic (Ap = 0) 
and insulating (Ap > 0) states. 

The Mott insulator in Fig. [l] exists for repulsive inter- 
actions {Uoo ^0)- It falls into the category COSl, with 
a nonzero charge gap Ap (and hence Agp > 0) reflecting 
the energy cost for doubly occupied sites, a vanishing spin 
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TABLE I: The different phases discussed in the text can be 
distinguished by the absence or presence of gaps for charge, 
spin, and single-particle excitations. The latter are denoted 
as Ap, Acj, and Asp, respectively. The last column gives the 
corresponding class in the Balents-Fisher notation. 



Phase 


A. 


A, 


Asp 


Class 


Luttinger liquid 


zero 


zero 


zero 


ClSl 


Luther-Emery liquid 


zero 


nonzero 


nonzero 


CISO 


Mott insulator 


nonzero 


zero 


nonzero 


COSl 


Peierls insulator 


nonzero 


nonzero 


nonzero 


COSO 



gap Act = 0, and dominant, power-law spin-density- wave 
(SDW) correlations. 

The intermediate metallic phase exists for small 
enough A and [/, with the explicit scale depending on 
^^^^I8i9iiiii7j one dimension, a metallic phase can ei- 
ther correspond to a Luttinger liquid with Agp = Ap = 
Act = (ClSl), or to a Luther-Emery liquid with Ap = 
but Asp, A^ > (ClSO).^The Luther-Emery liquid can 
be understood as a liquid of bosons, each consisting of 
two electrons bound into a spin singlet. The excitation 
spectrum for these bosons (corresponding to bipolarons 
in electron-phonon models) is gapless within the metal- 
lic phase (hence Ap = 0), but the binding of electrons 
into pairs gives rise to gaps for (electronic) single-particle 
and spin excitations (Agp , A^r > 0). The results pre- 
sented below suggest that in the Holstein and Holstein- 
Hubbard models, the metallic phase is a Luther-Emery 
liquid (CISO), except for the line Uoo = 0, which belongs 
to the Luttinger liquid class ClSl. 

The extent of the m etallic ph ase increases with increas- 
ing phonon f r eque ncy ) ^ I ^ I ^ ^ I ^ This can be ascribed to 
the suppression of charge order (the latter is generally 
favored for Uoo < 0) due to enhanced lattice fluctua- 
tions. These fluctuations have been shown^^ to destroy 
the Peierls state of the half-filled Holstein model below 
a critical value of A. The numerical results suggest that 
the metallic region exists in the Holstein-Hubbard model 
for Uoo < 0, and is adiabatically connected to that of 
the Holstein model {U = 0). For uJo/t = 5, the density- 
matrix renormalization group (DMRG)^ yields a larger 
metallic region than the QMC simulations,^ whereas for 
ooo/t = 0.5 the two methods are in quite good agreement. 
These two frequencies will also be considered here. Fi- 
nally, at the tricritical point where the Mott, Peierls and 
metallic phases intersect, a first-order phase transition 
has been reported.^ 

When the electron-phonon interaction dominates 
(^oo <^ 0), the system is a Peierls insulator with nonzero 
single-particle, charge and spin gaps (COSO). The origin 
of the spin gap is again a pairing of electrons into spin sin- 
glets, similar to the Luther-Emery phase. At T = 0, the 
Peierls state has long-range charge-density- wave (CDW) 
order with period q = 2A:f, corresponding to a pair of 
electrons residing at every other lattice site. As in the 
mean- field Peierls problem,'^^^ excitations out of the en- 



FIG. 2: (Color online) (a) Dynamical charge structure factor, 
(b) dynamical spin structure factor, and (c) single-particle 
spectral function for uJo/t = 5, U/t = 1 , and A = 0.5, cor- 
responding to the metallic phase (CISO).'^^ Here L — 30, 
/3t = 30. The dashed line indicates the charge velocity in the 
noninteracting case (A = 0, [/ = 0). 



ergetically favored 2/cf state cost a finite energy A^ 



A. Spectra in the metallic phase 

Figure [2] shows the dynamical charge and spin struc- 
ture factors, defined in Eqs. Q and d6|), as well as the 
single-particle spectral function [Eq. (|7f] in the nonadi- 
abatic regime uo/t = 5 and for U/t = 1^ X = 0.5. For 
these parameters, the metallic phase has been found to 
extend up to A ~ 0.65 by QMC simulations,^ and up to 
A ^ 0.9 by DMRG calculations. El 

The dynamical charge structure factor, plotted in 
Fig.|2|a), reveals a gapless mode at long wavelengths. Its 
velocity, Vp, is noticeably smaller than the noninteracting 
value Vp = 2t as a result of the enhanced mass of bipo- 
larons. At g = TT = 2/cf, there is an almost completely 
softened excitation with dominant spectral weight, which 
indicates strong but power-law charge correlations with 
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period 2/cf that are a precursor of the ordered Peierls 
phase. The dynamical spin structure factor in Fig. [2|b) 
also shows a linear mode. A possible spin gap is too small 
to be visible for these parameters. 

Figure [3] shows the excitation spectra for a stronger 
electron-phonon coupling A = 0.875, closer to (poten- 
tially on) the DMRG phase boundary for the metallic 
phase.^I^ The charge structure factor still appears gap- 
less, signaling metallic behavior, with an even smaller 
charge velocity and stronger softening ai q = 2/cf than 
in Fig. [2|a). The spin structure factor has developed a 
well visible gap at long wavelengths, which is also re- 
flected in the corresponding single-particle spectrum in 
Fig. [sj^c). The existence of bound singlets, as reflected 
by the spin gap, makes bipolarons the low-energy degrees 
of freedom. The corresponding Drude weight strongly 
depends on the phonon frequency.!^ An important point 
concerning Figs. [2] and [3] is that a gapless mode with sig- 
niflcant Drude weight exists in Sp{q^uo) near g = 0. In 
contrast, the spectrum in the insulating Peierls phase [see 
Fig. [7|^a)] exhibits a strong depletion of spectral weight 
and a flnite gap for long- wavelength charge excitations. 

The metallic region shrinks with decreasing phonon 
frequencyl ^ l -^-^ l -^^ l To study the adiabatic regime, a value 
uo^/t = 0.5 is considered, for which the extent of the 
metallic phase is known .^^^ For U /t = 0.2, it exists up to 
A ^ 0.25.^ Since the spin gap increases on approaching 
the Peierls phase (see discussion above for cjq/^ = 5), an 
electron-phonon coupling close to the phase boundary 
(A = 0.225) is chosen. 

The charge velocity Vp in Fig.|4|a) is almost unchanged 
compared to the noninteracting case, which can be under- 
stood in terms of the rather small bare coupling constant 
g. Hence, the electrons are only loosely bound into large 
bipolarons. An important difference to the case uo^/t = 5 
considered before is that the charge structure factor in 
Fig. [4ja) reveals the renormalized phonon frequency at 
low energies. The latter is partially softened at g = 2/cf, 
which is a previously noted precursor effect of the Peierls 
transition.aaa^ 

The spin structure factor in Fig. [4]^b) appears to be 
gapless (although it be will argued below that a flnite 
spin gap exists throughout the metallic phase), with a 
weakly renormalized, linear low-energy mode. Finally, 
the single-particle spectrum in Fig.[4]^c) also appears gap- 
less, with small but visible signatures of the hybridized 
polaron modes located near uo = ±uJo^^ 

Previous numerical calculations of excitation spectra 
of the half-fllled, metallic Holstein-Hubbard model were 
restricted to the single-particle spectrum, and showed 
in particular that spin-charge separation is only visible 
for large values of A (this in turn requires a large ratio 
oJo/t for a metallic state to exist), in accordance with sim- 
ulations at quarter fllling.^^ This conclusion explains the 
absence of such features in the present results. The spec- 
trum was also calculated at U = 0.'^ Analytical results 
include the exact spectrum of a Luttinger liquid (i.e., 
without a spin gap) coupled to phonons,^^ as well as ap- 



FIG. 3: (Color online) (a) Dynamical charge structure factor, 
(b) dynamical spin structure factor, and (c) single-particle 
spectral function for cjo/t = 5, U/t = 1, and A = 0.875, 
corresponding to the metallic phase (CISO)."'^"'^ Here L = 30, 
/3t = 20. The dashed line indicates the charge velocity in the 
noninteracting case (A = 0, [/ = 0). 



proximate results for the spectral function of the Luther- 
Emery liquid. The latter work predicts that branch cuts 
may be regularized by the existence of a gap. Interest- 
ingly, for small A where the spin gap is not resolved and 
where spin-charge separation is not detected in numer- 
ical simulations, the spectrum is very well described by 
analytical results for a spinless Luttinger liquid coupled 
to phonons,^^ and closely resembles the spectrum of the 
spinless Holstein model.^ 



B. Spin gap in the metallic phase 

A spin gap is visible in the dynamical correlation func- 
tions shown above, at least close to the Peierls phase 
boundary. The initially exponential dependence of this 
gap on the interaction strength and restrictions in system 
size make it diflicult to detect the spin gap at weak cou- 
pling by considering excitation spectra. Instead, evidence 
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FIG. 4: (Color online) (a) Dynamical charge structure factor, 
(b) dynamical spin structure factor, and (c) single-particle 
spectral function for uj^/t = 0.5, IJ jt — 0.2, and A = 0.225, 
corresponding to the metallic phase (CI SO). Here L = 50, 
/3t = 50. The dashed lines in (a) and (b) indicate the charge 
and spin velocities in the noninteracting case (A = 0, [7 = 0), 
respectively. The horizontal lines in (c) mark icjo- 



for a nonzero spin gap in the metallic Holstein-Hubbard 
model comes from numerical results for the Luttinger liq- 
uid parameter i^cr, which can in principle be determined 
from the spin structure factor and finite-size extrapola- 
tion as K(j = liiRL^oo ^S(j{qi) /qi, with qi = 27r/L. For 
the half-filled Hubbard model, were complications such 
as retarded electron-electron interactions are absent, nu- 
merical simulations generically give 7rSa{qi)/qi > 1 for 
U > 0, 7rSa{qi)/qi < 1 for < 0, and 7rSa{qi)/qi = 1 
exactly at = 0.^ In particular, although logarithmic 
corrections make it hard to demonstrate K^r = 1 or 
in the spin gapless and spin gapped phases, respectively, 
the finite-size estimates always decrease with increasing 
system size. Therefore, for the repulsive case U > 0, 
^Scr{qi) / qi slowly approaches the value 1 implied by spin 
rotation symmetry from above, whereas for the attrac- 
tive case /7 < 0, 7TScr{qi)/qi deviates more and more 
from 1. Within a low-energy theory, and given spin 
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FIG. 5: (Color online) Renormalized spin structure fac- 
tor TrSaiqi) /qi at qi = 27r/L, related to via — 
XmiL^ooT^Saiqi) /qi^ as a function of inverse system size in 
the metallic phase of the Holstein model. Results are for 
/3t = L, cjo/t = 0-5, = 0, and (a) A = 0.01, (b) A = 0.1. 



symmetry, the latter result can only be understood in 
terms of Ka^ = and a spin gap. Similar behavior can 
also be observed beyond the Hubbard model, and values 
^Scr{qi) / qi < 1 have been argued to represent empirical 
evidence for a spin gap.l^ For a more detailed discussion 
see Kefs. I9l58l 

Figure |5] shows the system size dependence of 
^Sa{qi) / qi in the adiabatic regime (cjqA = 0-^) of 
Holstein model in the metallic phase. Even for very weak 
coupling A = 0.01, 7TS(j{qi)/qi < 1 and monotonically 
decreasing with increasing system size. At stronger cou- 
pling A = 0.1, the size dependence is significantly more 
pronounced. A similar picture persists for larger phonon 
frequencies. The spin gap has also been measured di- 
rectly by means of the DMRG method. While A^r is 
clearly finite for selected parameters in the metallic re- 
gion, it is difficult to decide if this is true of the whole 
metallic phase. In particular, the exponential depen- 
dence of the gap on the coupling makes it practically 
impossible to detect the critical point directly from A^r- 

To provide further evidence for the existence of a spin 
gap, it is useful to consider the real- space correlation 
functions. In a gapless Luttinger liquid, the slowest de- 
caying correlations are given b}P^ 



^CDw(^) ~ cos(2/cFr)r , 
^SDw(0 - cos(2feFr)r-^^-^- , 
^SDw(0 - cos(2fepr)r-^^-i/^- 



(8) 



Rss{r) 



For a generic spin-rotation invariant system with K^r = 1, 
spin and charge correlations are degenerate, and domi- 
nate over pairing correlations if Kp < 1. On the other 
hand, attractive interactions {Kp < 1) lead to domi- 
nant pairing correlations. Additional complications are 
logarithmic corrections arising from marginally relevant 
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FIG. 6: (Color online) Real-space charge, pairing, and spin correlations (all rescaled by a factor 1/2) for (a) [/ = 0, A = 0.1, 
ujQ = 0.5t, (b) [7 = 0, A = 0.1, cjo = 5t, both corresponding to the metallic phase of the Holstein model, (c) The same 
correlation functions for the attractive Hubbard model with U — — 0.4t. (d) As in (a), but for U — ujq — 0.5t, and A = 0.2. 
Results are for /3t = L = 62. Spin correlations in the z and x directions are identical within error bars; S^^ is shown because 
of its smaller error bars. 



operators. In the repulsive Hubbard model, such cor- 
rections cause spin correlations to dominate over charge 
correlations.!^ As pointed out by Voit,^^ dominant q = 
2/cf CDW correlations cannot occur in a gapless Lut- 
tinger liquid. (Although logarithmic corrections could in 
principle favor charge over spin correlations, there is no 
known example of such behavior. Moreover, such correc- 
tions would not explain K(j < 1.) 

In contrast, due to the exponential suppression of spin 
correlations, a Luttinger liquid with a gap for spin ex- 
citations supports dominant 2/cf charge correlations. In 
terms of the low-energy model, this case arises when at- 
tractive backscattering is taken into account. The cor- 
responding model is often referred to as a Luther-Emery 
model. A familiar example is the attractive Hubbard 
model. That attractive backscattering can originate from 
electron-phonon interaction is illustrated by the previ- 
ously mentioned explicit relation between the Holstein 
model and the attractive Hubbard model. To obtain 
the correlation functions of the a Luther-Emery liquid, 
the spin gap can formally be accounted for by setting 
Kcj = 0, leading tcF^ i?CDw(^) ~ cos(2/cf^)^~^^ and 
Rss{r) ~ ^-i/^p^ also Eq. Hence, neglecting 

possible logarithmic corrections, repulsive (attractive) in- 
teractions lead to dominant charge (pairing) correlations. 

Results for the charge, spin and pairing correlation 
functions are shown in Fig. [gJ Figures [6]^a) and (b) cor- 



respond to the adiabatic {uJo/t = 0.5) and the nonadi- 
abatic regime {uo/t = 5), respectively. To be able to 
explore very weak interactions Uoo, it is advantageous to 
consider U = 0, so that the metallic phase exists down 
to A = 0. This choice also eliminates any uncertainties 
about the phase boundaries of the metallic region. For 
the coupling A = 0.1 chosen, charge correlations domi- 
nate over pairing and spin correlations, as expected for 
a Luther-Emery liquid with repulsive interactions. This 
dominance is more pronounced for uj^/t = 0.5 [Fig. [6]^ a)] 
than for uo/t = 5 [Fig. [6]^b)] because the system is closer 
to the Peierls phase, and cannot be explained by a rescal- 
ing of the constant (i.e., independent of r) prefactors of 
2/cf correlations in the charge and spin channels. The re- 
sults are qualitatively similar also at weaker interactions 
(the values investigated where as small as A = 0.01), but 
it becomes increasingly more difficult to distinguish the 
different correlation functions on finite systems. 

For comparison. Fig. ^c) presents the same correla- 
tion functions for the attractive Hubbard model with 
U — —XW — — 0.4t, corresponding to the limit uq = 
oo at the same coupling A = 0.1 used in Figs. [6](a), 
(b). Whereas spin correlations are again clearly sup- 
pressed, the mapping between the attractive and the re- 
pulsive Hubbard model at half filling implies a degener- 
acy of charge and pairing correlations that is captured, 
within statistical errors, by the numerical data. Finally, 
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FIG. 7: (Color online) (a) Dynamical charge structure factor, 
(b) dynamical spin structure factor, and (c) single-particle 
spectral function for uJo/t = 0.5, IJ jt — 0.2, and A = 0.4, 
corresponding to the Peierls insulator (COSO). Here L = 30, 
/3t = 30. 



Fig. ^d) shows results at a stronger coupling A = 0.2, 
for which charge correlations and the suppression of spin 
correlations are much more pronounced. The exponen- 
tial decay of spin correlations has been observed before 
in DMRG calculations.^^ 

The results for 7tScj{qi)/qi and the real-space correla- 
tion functions are fully compatible with the existence of 
a spin gap throughout the metallic phase. In particular, 
the behavior of these observables is essentially identical 
to the attractive Hubbard model, for which the existence 
of a spin gap is well established. In contrast to the gapless 
Luttinger liquid model, the Luther-Emery model hence 
provides a consistent description of the numerical data. 
The existence of a spin gap for any nonzero electron- 
phonon coupling in the Holstein model has also been pre- 
dicted based on renormalization group calculations .^^^^ 

The above results suggest that the half-filled, metal- 
lic Holstein model behaves as a repulsive Luther-Emery 
liquid, with dominant charge correlations, and a nonzero 
spin gap in the thermodynamic limit that suppresses spin 



correlations with respect to charge and pairing correla- 
tions. Assuming that the metallic phase is the same 
at = and U > leads to the conclusion that the 
whole metallic region of the Holstein-Hubbard model is 
a Luther-Emery liquid. The exponential length scale re- 
lated to the spin gap is expected to lead to a crossover 
(as a function of system size) from the correlation func- 
tions ([8|, as appropriate for a system with gapless spin 
excitations, to the corresponding analytical results for 
the Luther-Emery liquid.— In particular, the relation 
<^CDW<^ss = I7 expected for a Luther-Emery liquid, is not 
obeyed by the numerically determined power-law expo- 
nents. Whereas this fact has previously been attributed 
to retardation effects,^ it can be expected to be a general 
problem associated with finite-size simulations of Luther- 
Emery liquids. 



C. Spectra in the Peierls phase 

The single-particle spectral function in the Peierls 
phase h as been calculated before using numerical 
methods.f^^ Bosonization results for the spectrum of a 
quarter-filled CDW insulator are also available. The dy- 
namical charge structure factor has been studied for the 
spinless Holstein model.S^ Using the CTQMC method, 
the Peierls state is most accessible in the adiabatic 
regime. Specifically, ujo/t = 0.5, A = 0.4, and U/t = 0.2 
are considered, and results are shown in Fig. ^ 

The dynamical charge structure factor in Fig. [7|^a) is 
characterized by a finite gap and strongly suppressed 
spectral weight for long- wavelength charge excitations, 
and by a soft phonon mode with gapless excitations at 
q = 2/cf. Similar to the metallic Luther-Emery phase, 
the Peierls state has a finite (but larger) spin gap, see 
Fig. [rj^b). Compared to the spin structure factor in the 
metallic phase. Fig. |4jb), there is also a finite gap for 
spin excitations with q = 2/cf- 

The single-particle spectrum in Fig. [tJ^c) reveals the 
expected Peierls gap, as well as clear signatures of back- 
folded shadow bands at high energies — as a result of 
dimerizatioiP^I — and soliton excitations at low energies. 
The latter have previously been observed in the spin- 
less Holstein model, and the extended Holstein model 
with nonlocal interactions,^ and distinguish the otherwise 
qualitatively (apart from the size of the single-particle 
gap) similar spectra of the metallic phase [Fig. ^c)] and 
the Peierls phase [Fig. [tJ^c)]. 



D. Spectra in the Mott phase 

The excitation spectra in the Mott phase are shown 
in Fig. [sj for the parameters ujo/t — 5, A = 0.25, and 
U/t = A. The single-parti cle spect rum in this phase has 
previously been calculated. ^^^^^SIM] 

In accordance with the discussion of the phase dia- 
gram at the beginning of Sec. |IV[ the results in Fig. [S] 
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FIG. 8: (Color online) (a) Dynamical charge structure factor, 
(b) dynamical spin structure factor, and (c) single-particle 
spectral function for uj^/t = 5, U/t = 4, and A = 0.25, corre- 
sponding to the Mott insulator (COSl). Here L = 30, /3t = 20. 



reveal gapped charge excitations but gapless spin exci- 
tations, i.e., just opposite to the metallic phase shown 
in Fig. [3j Although the spin symmetry is not broken in 
the Mott phase, the spin structure factor reveals gapless 
spin excitations ai q = 2/cf reminiscent of the soft phonon 
mode related to long-range charge order in Fig.[7|^a). For 
the value uJo/t = 5 considered, the renormalized phonon 
mode is not visible in Sp{q,uj), and it is strongly sup- 
pressed compared to the Peierls state even for ujo/t = 0.5 
(data not shown). 

The single-particle spectrum in Fig. [sj^c) bears a close 
resemblance to the results for the Hubbard model with 
the same U/t = 4,^^ although signatures of spin-charge 
separation are suppressed here as a result of a reduced, 
effective interaction {Uoo/t = 3) and thermal fluctuations 
at T > 0. The effects of the latter are quite subtle for 
intermediate values of U/t.^^ Instead of the pronounced 
backfolded shadow bands related to long-range charge or- 
der in the Peierls state, Fig.[7|^c), there is rather incoher- 
ent spectral weight with no clear dispersion. Addition- 



ally, the hallmark soliton excitations visible in Fig. [tJ^c) 
are completely absent. 



V. CONCLUSIONS 

The single- and two-particle excitation spectra of 
the one-dimensional, half-filled Holstein-Hubbard model 
have been calculated with the continuous-time quantum 
Monte Carlo method. The spectra in the metallic phase 
are consistent with a Luther-Emery liquid that has gap- 
less charge excitations but gapped spin excitations with 
a gap that opens exponentially as a function of the inter- 
action strength. In the Peierls phase, the spectra reveal 
both a charge and a spin gap, and a soft charge mode 
dit q = 2/cf. The single-particle spectral function reveals 
backfolded shadow bands and soliton excitations. The 
Mott phase has a charge gap but no spin gap, as well 
as a soft spin mode despite the absence of long-range or- 
der. No clear signatures of spin-charge separation were 
observed. 

The static charge and spin correlation functions have 
been calculated on rather large systems and reveal dom- 
inant 2/cf charge correlations even for very weak cou- 
pling. Because such behavior cannot occur in a generic 
Luttinger liquid, this observation can be regarded as evi- 
dence for a nonzero spin gap.^^ In particular, the behav- 
ior of these correlation functions is essentially the same 
as for the attractive Hubbard model (for which the ex- 
istence of a spin gap is generally accepted) even for low 
phonon frequencies. When combined with the fact that 
the Luttinger parameter K^j is found to be less than unity 
even for very weak coupling, the numerical data are natu- 
rally explained by assuming that the metallic phase of the 
Holstein-Hubbard model is a Luther-Emery liquid, with 
electrons paired into singlet bipolarons. E arlier reports of 
both gapless and gapped metallic region^HEH c^n be at- 
tributed to the exponentially small size of the gap which 
makes its detection very demanding. The existence of 
a nonzero spin gap also has important implications for 
the low-energy description of this model. In particular, 
the spin gap is not captured by the exact results avail- 
able for the spectral function of the Luttinger model with 
phonons.^^ Finally, given the relevance of backscattering 
at any band filling, the Luther-Emery description can be 
expected to be generic for metallic phases in Holstein- 
type models, in accordance with the suggestive physical 
picture of electron pairs that are bound into bipolarons. 



Acknowledgments 

Computer time at the Jiilich Supercomputing Cen- 
tre and financial support from the DEC (Grant 
No. Ho 4489/2-1) are gratefully acknowledged, as are dis- 
cussions with S. Ejima, H. Fehske, V. Meden, A. Sandvik, 
D. Schuricht, S. Wessel, and, in particular, F. Essler. 



10 



^ B. Sipos, A. F. Kusmartseva, A. Akrap, H. Berger, 
L. Forro, and E. Tutis, Nature Mat. 7, 960 (2008). 

^ M. Hohenadler, F. F. Assaad, and H. Fehske, Phys. Rev. 
Lett. 109, 116407 (2012). 

^ M. Hohenadler and W. von der Linden, Phys. Rev. B 71, 
184309 (2005). 

^ J. T. Devreese and A. S. Alexandrov, Rep. Prog. Phys. 72, 
066501 (2009). 

^ F. Guinea, J. Phys. C: Sohd State Phys. 16, 4405 (1983). 
^ J. E. Hirsch, Phys. Rev. B 31, 6022 (1985). 
^ J. Stein, Europhys. Lett. 39, 413 (1997). 
^ R. T. Clay and R. P. Hardikar, Phys. Rev. Lett. 95, 096401 
(2005). 

^ R. P. Hardikar and R. T. Clay, Phys. Rev. B 75, 245103 
(2007). 

^° H. Fehske, G. Wellein, G. Hager, A. WeiBe, and A. R. 

Bishop, Phys. Rev. B 69, 165115 (2004). 

H. Fehske, G. Hager, and E. Jeckelmann, Europhys. Lett. 

84, 57001 (2008). 
^2 S. Ejima and H. Fehske, J. Phys. Conf. Ser. 200, 012031 

(2010) . 

W. Q. Ning, H. Zhao, C. Q. Wu, and H. Q. Lin, Phys. Rev. 
Lett. 96, 156402 (2006). 

H. Matsueda, T. Tohyama, and S. Maekawa, Phys. Rev. B 
74, 241103 (2006). 

K.-M. Tarn, S.-W. Tsai, and D. K. Campbell, Phys. Rev. 
B 84, 165123 (2011). 

M. Tezuka, R. Arita, and H. Aoki, Phys. Rev. Lett. 95, 
226401 (2005). 

M. Tezuka, R. Arita, and H. Aoki, Phys. Rev. B 76, 155114 
(2007). 

A. Payeur and D. Senechal, Phys. Rev. B 83, 033104 

(2011) . 

S. Reja, S. Yarlagadda, and P. B. Littlewood, Phys. Rev. 
B 84, 085127 (2011). 

S. Reja, S. Yarlagadda, and P. B. Littlewood, Phys. Rev. 

B 95, 096401 (2012). 
2^ T. Holstein, Ann. Phys. (N.Y.) 8, 325; 8, 343 (1959). 
'^'^ J. E. Hirsch and E. Fradkin, Phys. Rev. B 27, 4302 (1983). 

E. Jeckelmann, C. Zhang, and S. R. White, Phys. Rev. B 

60, 7950 (1999). 

L P. Bindloss, Phys. Rev. B 71, 205113 (2005). 
2^ K.-M. Tam, S.-W. Tsai, D. K. Campbeh, and A. H. Castro 
Neto, Phys. Rev. B 75, 161103 (2007). 
D. Meyer, A. C. Hewson, and R. Bulla, Phys. Rev. Lett. 
89, 196401 (2002). 

W. Roller, D. Meyer, Y. Ono, and A. Hewson, Euro. Phys. 
Lett. 66, 559 (2004). 
2^ W. Roller, D. Meyer, and A. C. Hewson, Phys. Rev. B 70, 
155103 (2004). 

M. Capone, G. Sangiovanni, C. Castellani, C. Di Castro, 

and M. GrilU, Phys. Rev. Lett. 92, 106401 (2004). 

P. Werner and A. J. Millis, Phys. Rev. Lett. 99, 146404 

(2007). 

J. Bauer and A. C. Hewson, Phys. Rev. B 81, 235113 
(2010). 

R. Peierls, Surprises in Theoretical Physics (Princeton 
University Press, New Jersey, 1979). 

G. Beni, P. Pincus, and J. Ranamori, Phys. Rev. B 10, 



1896 (1974). 

H. Bakrim and C. Bourbonnais, Phys. Rev. B 76 (2007). 
Y. Takada and A. Chatterjee, Phys. Rev. B 67, 081102 
(2003). 

H. Fehske, A. P. Rampf, M. Sekania, and G. Wellein, Eur. 
Phys. J. B 31, 11 (2003). 

M. Hohenadler and F. F. Assaad, J. Phys.: Condens. Mat- 
ter 25, 014005 (2013). 

A. Luther and V. J. Emery, Phys. Rev. Lett. 33, 589 
(1974). 

V. J. Emery, in Highly conducting one- dimensional solids^ 
edited by J. T. Devreese (Plenum Press, New York, 1979), 
p. 247. 

J. Voit, Eur. Phys. J. B 5, 505 (1998). 
T. Giamarchi, Quantum physics in one dimension (Claren- 
don Press, Oxford, 2004). 

R. Citro, E. Orignac, and T. Giamarchi, Phys. Rev. B 72, 
024434 (2005). 

M. Hohenadler, D. Neuber, W. von der Linden, G. Wellein, 

J. Loos, and H. Fehske, Phys. Rev. B 71, 245111 (2005). 

M. Hohenadler, G. Hager, G. Wellein, and H. Fehske, J. 

Phys.: Condens. Matter 19, 255202 (2007). 

J. Voit, Rep. Prog. Phy. 58, 977 (1995). 

F. F. Assaad, Phys. Rev. B 78, 155124 (2008). 

H. Fehske, G. Wellein, A. Weisse, F. Gohmann, H. Buttner, 

and A. R. Bishop, Physica B 312-313, 562 (2002). 

A. N. Rubtsov, V. V. Savkin, and A. L Lichtenstein, Phys. 

Rev. B 72, 035122 (2005). 

M. Hohenadler, H. Fehske, and F. F. Assaad, Phys. Rev. 
B 83, 115105 (2011). 
^° F. F. Assaad and T. C. Lang, Phys. Rev. B 76, 035116 
(2007). 

E. Guh, A. J. Millis, A. L Lichtenstein, A. N. Rubtsov, 
M. Troyer, and P. Werner, Rev. Mod. Phys. 83, 349 (2011). 

^2 R. S. D. Beach, arXiv:0403055 (2004). 

L. Balents and M. P. A. Fisher, Phys. Rev. B 53, 12133 
(1996). 

M. Hohenadler, G. Wellein, A. R. Bishop, A. Alvermann, 
and H. Fehske, Phys. Rev. B 73, 245120 (2006). 
C. E. Creffield, G. Sangiovanni, and M. Capone, Eur. Phys. 
J. B 44, 175 (2005). 

S. Sykora, A. Hiibsch, R. W. Becker, G. Wellein, and 
H. Fehske, Phys. Rev. B 71, 045112 (2005). 
V. Meden, R. Schonhammer, and O. Gunnarsson, Phys. 
Rev. B 50, 11179 (1994). 

R. T. Clay, A. W. Sandvik, and D. R. Campbeh, Phys. 
Rev. B 59, 4665 (1999). 

H. J. Schulz, Phys. Rev. Lett 64, 2831 (1990). 

R. Yonemitsu and M. Imada, Phys. Rev. B 54, 2410 

(1996). 

F. H. L. Essler and A. M. Tsvelik, Phys. Rev. Lett. 88, 
096403 (2002). 

J. Voit, L. Perfetti, F. Zwick, H. Berger, G. Margaritondo, 

G. Griiner, H. Hochst, and M. Grioni, Science 290, 501 
(2000). 

®^ D. Senechal, D. Perez, and D. Plouffe, Phys. Rev. B 66, 
075129 (2002). 



